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Abstract 

\, 

Seemingly  identical  Nd:YAG  laser  rods  often  exhibit  sig- 
nificant variations  in  output  power.  The  purpose  of  this  re- 
port is  to  investigate  the  possible  cause  or  causes  for  the 
performance  variations  in  Nd:YAG  laser  rods. 

The  fundamental  chemical  and  crystal  growth  characteris- 
tics of  Nd:YAG  are  reviewed.  The  spectroscopy  of  Nd:YAG  is 
discussed  with  emphasis  on  the  laser  transition.  The  more 
common  spectroscopic  parameters  are  defined  and  the  important 
relationships  between  them  derived.  Expressions  relating  the 
laser  performance  of  CW  and  Q-switched  systems  to  the  host 

parameters  are  developed  and  the  important  parameters 

? 

jPg-pT— are  discussed.  Possible  energy  transfer  mechanisms 
to  account  for  performance  variations  are  presented  and  eval- 
uated. 

It  is  proposed  that  concentration  variations  are  primar- 
ily responsible  for  laser  performance  fluctuations  and  a meth- 
od for  obtaining  experimental  verification  is  outlined. 


Vi 


OUTPUT  POWER  VARIATIONS 


IN 

NDrYAG  LASER  HOSTS 


I.  Introduction 


In  1958,  Schawlow  and  Townes  first  proposed  the  laser 
(Ref  1).  In  i960,  Maiman  demonstrated  the  first  successful 
laser  system  employing  a ruby  crystal  in  the  pulsed  mode  of 
operation  (Ref  2).  Following  the  success  of  Maiman,  laser 
action  was  achieved  using  a wide  variety  of  active  ions  and 
host  materials. 

In  1961,  Johnson  and  Nassau  (Ref  3)  obtained  laser  ac- 
tion from  neodymium  (Nd)  In  calcium  tungstate  (CaWO^)  and 
Snltzer  (Ref  4)  observed  lasing  from  neodymium-doped  glass. 
Over  the  next  few  years,  neodymium  was  lased  in  numerous  host 
materials.  Then,  In  1964,  Gueslc,  et  al.,  (Ref  5)  demonstra- 
ted laser  action  In  neodyralum-doped  yttrium  aluminum  garnet 
(Nd:YAG).  In  the  years  since  Its  initial  discovery,  Nd:YAG 
has  emerged  as  one  of  the  more  popular  solid  state  laser  mat- 
erials (Ref  6:44). 

However,  Nd:YAG  Is  not  without  Its  drawbacks.  When 
operated  at  high  continuous  wave  (CW)  powers,  cherraally  In- 
duced focusing  and  birefringence  effects  produce  decreased 
polarized  and  TEMqq  power  outputs  (Ref  7,8).  For  neodymium 
concentrations  In  excess  of  1.2  atom  percent,  the  fluorescent 
lifetime  decreases  rapidly  resulting  In  Increased  threshold 
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pumping  requirements  CRef  9;  205).  Tlie  above  problems  repre- 
sent limitations  on  the  operation  of  Nd:YAG  laser  systems. 

I 

However,  even  if  a NdiYAG  system  is  operated  within  these  lim- 
itations, another  problem  exists;  seemingly  identical  laser 
rods  will  often  exhibit  significant  variations  in  power  out- 
put. The  cause (s)  of  these  power  variations  is  not  known. 

Thus,  the  primary  purposes  of  this  investigation  are  to 
propose  possible  mechanisms  for  the  observed  power  variations 
among  Nd:YAG  laser  rods  and  future  areas  of  investigation  to 
aid  in  determining  the  mechanism(s)  responsible  for  those 
variations.  Additionally,  it  is  hoped  that  this  thesis  will 
serve  as  an  Introductory  source  document  for  those  Interested 
in  Nd:YAG  materials  research. 

In  order  to  accomplish  the  above,  the  chemistry  and  crys- 
tal growth  of  Nd:YAG  are  discussed  briefly  in  Chapter  II. 

Then,  in  Chapter  III,  the  spectroscopy  of  Nd:YAG  is  discussed 
along  with  the  fundamental  spectroscopic  parameters  and  how 
they  are  determined.  In  Chapter  IV,  the  Important  equations 
for  the  four  level  laser  are  developed  and  the  key  host  para- 
meters identified.  Possible  mechanisms  for  power  variations 
are  presented  in  Chapter  V.  Recommended  areas  for  future  in- 
vestigations are  discussed  in  Chapter  VI.  A supplementary 
bibliography,  for  those  Interested  in  more  detailed  discussions 
of  the  topics  covered  in  this  thesis,  is  presented  in  Appen- 
dix A. 
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II.  The  Chemistry  and  Crystal 
Growth  of  Nd:YAG 


The  Chemistry  of  NdiYAG 

Neodymium  Is  a part  of  the  lanthanide  series  of  elements 
which  are  more  commonly  referred  to  as  the  rare-earths . The 
rare-earth  atoms  have  electronic  configurations  given  by: 

Is^  2s^  2p®  3s^  3p®  ‘’^s^4p®  4d^  “ 4f'^5s^5p®  5d”'6s^ 

where  m=0,l  and  n=0 ,2 ,3 9 • • . >1^ . They  are  primarily  trlvalent , 
although  some  of  the  atoms  may  also  be  divalent  or  tetravalent 
(Ref  10:15).  In  the  trlvalent  form,  the  electronic  config- 
uration Is  that  shown  above  except  that  the  5<i®  and  6s ^ elec- 
trons are  missing.  Thus,  the  outermost  electrons  (5s^5p®) 
comprise  the  xenon  rare-gas  shell  and  are  optically  Inactive. 
Inside  the  xenon  shell  lies  the  4f'^  shell  with  Its  n optically 
active  electrons  (n=0,l, . . . ,14) . Since  the  4f  shell  lies  In- 
side the  xenon  shell.  Its  electrons  Interact  weakly  with  other 
Ions  resulting  In  sharp  spectral  lines,  high  fluorescent 
yields,  and  remarkably  similar  absorption  and  emission  spectra 
for  a given  Ion  In  various  hosts  (Ref  11:267).  The  neodymium 
atom  (atomic  number=60,  atomic  welght=l44 . 24 ) corresponds  to 
the  m*0,  n=4  electronic  configuration,  while,  for  the  trlva- 
lent neodymium  Ion,  n=3. 

Yttrium  aluminum  garnet  has  the  molecular  formula, 

Y3Al50i2>  with  molecular  weight,  593.7.  The  overall  crystal 
structure  exhibits  cubic  symmetry  with  eight  molecules  per 
unit  cell  and  a cube  edge  dimension  of  12  angstroms  (Ref  12:87). 
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Because  of  Its  cubic  symmetry,  YAG  exhibits  isotropic  thermal 
and  optical  behavior.  YAG  has  a high  coefficient  of  thermal 
conductivity  and  a low  thermal  expansion  coefficient,  both  of 
which  are  highly  desirable  in  a laser  host.  Additionally, 

YAG  is  very  hard  which  allows  it  to  be  machined,  polished 
and  handled  without  fear  of  fracture  (Ref  9:202). 

When  YAG  is  doped  with  neodymium,  the  neodymium  ions  re- 
place the  yttrium  ions  in  the  garnet  structure  resulting  in  a 
molecular  formula  of  Yj_^Nd^Al ^0^ ^ j where  x represents  the 
number  of  yttrium  ions  replaced  by  neodymium.  For  example, 
if  the  neodymium  concentration  is  one  atom  percent,  x=0.03 
(Ref  13).  For  normal  neodymium  concentrations  (1. 0-1.2  atom 
percent),  the  properties  of  Nd:YAG  are  essentially  unchanged 
from  those  of  pure  YAG.  Excellent  summaries  of  the  important 
properties  of  YAG  and  Nd:YAG  may  be  found  in  Ref  6 and  Ref 
12:88-92. 

The  Crystal  Growth  of  Nd : YAG 

Although  Nd:YAG  crystals  have  been  grown  by  several  dif- 
ferent methods,  the  preferred  growth  technique  is  Czochralski 
pulling.  Crystals  are  grown  by  this  method  by  dipping  a rota- 
ting seed  crystal  into  a crucible  of  molten  material  and  with- 
drawing it  at  a constant  rate.  The  growth  of  Nd:YAG  requires 
the  use  of  an  iridium  crucible  due  to  the  high  melting  point 
(1970  °C)  of  Nd:YAG.  Use  of  the  iridium  crucible  requires  that 
the  growth  process  be  conducted  in  an  inert  atmosphere  (typi- 
cally argon  and  nitrogen)  to  avoid  oxidation  of  the  crucible 
and  subsequent  contamination  of  the  melt  (Ref  12:85).  The 


seed  crystal  Is  normally  rotated  at  40-50  rpm  and  withdrawal 
rates  usually  range  from  0.1  to  1.5  mm/hr  depending  on  the 
dopant  level  and  crystal  (or  boule)  diameter.  The  growth  of 
high  quality  Nd:YAG  requires  that  the  starting  oxides  (YjOj, 
AljOj,  NdjOj)  used  in  the  melt  be  dry  and  of  99-999  percent 
purity  (Ref  9*. 203)-  It  Is  also  Important  that  the  boule  dia- 
meter and  the  shape  of  the  boule-melt  Interface  be  rigidly 
controlled  (Ref  l4:269). 

Even  under  optimum  growth  conditions,  Nd:YAG  boules  exhi- 
bit some  undesirable  growth  characteristics.  Three  of  these, 
core  effect,  strain,  and  melt  enrichment  are  important  In  the 
production  of  laser  rods  and  are  discussed  below. 

Core  Effect . During  Czochralskl  growth,  facets  form  at 
the  crystal  Interface  along  either  the  {211}  or  {110}  planes 
depending  on  the  growth  direction.  The  net  effect  of  these 
facets  Is  to  produce  an  optically  inhomogeneous  core  that  runs 
the  length  of  the  boule  and  Is  unsuitable  for  laser  rods. 
Although  this  phenomena  Is  well  known,  the  exact  cause  Is  not. 
One  proposal  Is  that  It  Is  due  to  a difference  In  the  Impurity 
segregation  on  the  planar  faceted  and  non-faceted  Interfaces 
(Ref  9:203).  Vhlle  the  core  diameter  (1-2  mm)  appears  to  be 
Independent  of  the  boule  diameter  (Ref  14:269),  It  does  vary 
with  growth  direction  with  the  smallest  core  resulting  from 
growth  In  the  <111>  direction  (Ref  9:204). 

Strain.  Boules  grown  In  the  <111>  direction  have  a 
roughly  hexagonal  cross  section  due  to  the  faceting  of  the 
{211}  planes.  When  the  boule  Is  viewed  In  a Twyman-Green 
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interferometer,  areas  of  high  strain  corresponding  to  the  facet 
boundaries  are  quite  evident  as  illustrated  In  Fig.  1.  It  Is 


Important  that  laser  rods  be  cut  only  from  the  strain  free 
areas  of  the  boule. 

Melt  Enrichment . Since  the  neodymium  'ion  Is  larger  than 
the  yttrium  Ion  it  replaces,  neodymium  tends  to  be  retained 
in  the  melt.  A measure  of  this  tendency  Is  given  by  the  ini- 
tial distribution  coefficient,  kj^,  defined  by 

ki  . C^/C^  (1) 

where  is  the  Initial  Nd  concentration  in  the  crystal,  and 
Is  the  initial  Nd  concentration  in  the  melt.  For  NdrYAG, 
the  value  of  ranges  from  0.12  (Ref  14:269)  to  0.21  (Ref  15). 
Thus,  as  the  boule  is  pulled  from  the  melt,  the  concentration 
of  neodymium  in  the  melt  Increases  and  the  melt  is  said  to  be 
enriched.  The  Increased  neodymium  concentration  in  the  melt 
leads  to  an  Increase  in  the  neodymium  concentration  in  the 
boule  as  it  is  pulled  from  the  melt.  The  resulting  concentra- 
tion gradient  can  be  predicted  by  the  normal  freezing  equation; 
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(2) 


<=l  ’ C,k^(l-g)''l-l 

where  * Nd  concentration  in  the  boule  at  length  *, 

® initial  Nd  concentration  in  the  melt 
* initial  distribution  coefficient 
g « fraction  of  the  melt  pulled  corresponding  to 
length  1. 

Figure  2 shows  the  variation  of  as  a function  of  g. 


Pig.  2.  Impurity  Concentration  Profile 
for  Normal  Freezing  (after  Ref  16:222) 


Experiments  by  Belt,  et  al.,  have  confirmed  the  validity  of 
equation  (2)  for  Nd:YAG.  Increases  of  20-30  percent  in  neo- 
dymium concentration  along  the  length  of  a 20  cm  boule  are 
typical.  Thus,  for  one  atom  percent  initial  neodymium  concen- 
tration, end  to  end  variations  along  a 30-50  mm  rod  are  on  the 
order  of  0.05-0.10  atom  percent  and  the  maximum  rod  to  rod 
variations  among  rods  from  the  same  boule  are  0.10-0.15  atom 
percent  (Ref  14:269-270). 

Melt  enrichment  places  serious  limitations  on  the  amount 
of  material  that  can  be  crystallized  from  the  melt.  For  exam- 
ple, if  the  initial  neodymium  concentration  in  the  boule  is 
one  atom  percent,  when  20  percent  of  the  melt  has  been  crys- 
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tallzed  the  concentration.  In  the  boule  has  Increased  to  ap- 
proximately 1.2  atom  percent  - the  threshold  for  concentration 
quenching.  When  40  percent  of  the  melt  has  been  crystallzed, 
the  boule  concentration  has  Increased  to  roughly  1.5  atom 
percent  and  the  growth  of  high  optical  quality  Nd:YAG  becomes 
extremely  difficult  (Ref  9:205).  Thus,  for  the  more  common 
doping  levels  only  30-40  percent  of  the  melt  can  be  crystal- 
lzed. 
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III.  Spectroscopy  of  Nd ; YAG 


Knowledge  of  emission  and  absorption  spectra  is  essen- 
tial in  the  study  of  any  laser  material.  In  this  section, 
the  spectra  for  Nd:YAG  are  discussed  with  emphasis  on  those 
areas  relating  to  the  performance  of  Nd:YAG  as  a laser  mater- 
ial. Also,  the  fundamental  spectroscopic  parameters  are  de- 
fined and  methods  of  determining  them  discussed.  Finally, 
various  values  of  selected  spectroscopic  parameters  for  NdiYAG 
are  presented  and  the  discrepancies  among  them  discussed. 

Absorption  and  Emission  in  Nd; YAG 

The  optical  and  near-infrared  absorption  and  emission 
observed  in  Nd:YAG  are  primarily  due  to  transitions  between 
energy  states  within  the  4f  subshell.  Now,  4f->4f  electric 
dipole  transitions  are  not  allowed  in  the  free  ion  by  Laporte's 
rule  since  the  4f  wavefunctions  all  have  the  same  parity. 
However,  placing  neodymium  ions  in  the  crystal  field  of  YAG 
causes  the  4f  wavefunctions  to  have  mixed  parities  leading 
to  the  so-called  forced  electric  dipole  or  Laporte  forbidden 
transitions  between  the  4f  states  (Ref  10:118). 

The  energy  level  diagram  for  Nd:YAG  can  be  derived  from 
the  absorption  and  emission  spectra.  A partial  energy  level 
diagram  for  NdiYAG  is  presented  in  Fig.  3a  with  the  energy 
level  term  designators  based  on  the  Russell-Saunders  coupling 
scheme.  The  most  intense  absorption  bands  Involve  transitions 
from  the  ground  state  C’lj/j)  to  the  "'P s/2'^^^9/2  ^ind  j/2 

levels  which  correspond  to  wavelength  ranges  of  7300-7700  ang- 


Pig.  3.  Energy  Levels  in  Nd:YAG.  (a)  Partial  Energy 
Level  Diagram  for  Nd  In  YAG  at  4.2  (Ref  19:A712). 
(b)  Stark  Splittings  of  the  ‘‘P,/,  and  ‘*1  manifolds  at 
77  °K  (Ref  20:451). 
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Stroms  and  7900-8300  angstroms  CRef  17:3^-35).  According  to 
Brandewie  and  Talk  (Ref  l8),  all  the  levels  below  the  “*^3/2 
(3550  angstroms)  relax  nonradlatively  to  the  upper  laser 
level  C’Ps/z),  therefore,  the  best  pump  sources  for  Nd:YAG 
are  those  with  large  fractions  of  their  spectral  output  In 
the  7300-8300  angstrom  range. 

Since,  for  the  typical  excitation  wavelengths  used  with 
Nd:YAG,  all  the  fluorescent  emissions  originate  from  the  “Rj/j 
manifold,  the  relevant  energy  levels  for  the  emission  spec- 
trum are  those  shown  In  Pig.  3b.  As  Indicated  In  that  figure, 
the  (2J+l)-fold  degeneracy  of  the  various  manifolds  predicted 
by  the  Russell-Saunders  coupling  for  the  free  Ion  has  been 
partially  removed  due  to  the  effects  of  the  crystal  field. 

The  resulting  levels  are  referred  to  as  Stark  levels  and.  In 
Pig.  3b,  are  labeled  according  to  the  convention  used  by  Dleke 
(Ref  10).  Por  Nd:YAG,  each  of  the  Stark  levels  Is  Itself 
doubly  degenerate  (Ref  21:2569).  Pluorescent  emission  from 
the  **1  manifolds  has  not  been  observed  and  It  Is  assumed  that 
these  levels  rapidly  relax  nonradlatively  to  the  ground  state. 
As  Indicated  by  the  arrows  In  Pig.  3b,  the  laser  transition 
(X*1.06u)  for  Nd:YAG  Is  actually  a composite  of  the  Rj'^j  and 
Rj-Yj  transitions  between  the  '*^3/2  '’^11/2  nianlfolds  (Ref 

22). 

Pundamental  Spectroscopic  Parameters 

In  the  study  of  laser  materials,  certain  spectroscopic 
parameters  are  frequently  used  to  describe  the  merits  of  a 
material  as  a laser  medium.  The  more  Important  of  these  para- 
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meters  are  discussed  below. 


Llneshape.  The  dominant  line-broadening  mechanism  in 
Nd:YAG  is  thermal  broadening  (a  homogeneous  mechanism),  there- 
fore, the  llneshape  for  Nd:YAG  Is  Lorentzlan  (Ref  23:362). 

The  nonnallzed  Lorentzlan  llneshape  function  can  be  expressed 
as  : 

Av 

g(v)  = — {(v-v  )^  + (Av/2)*}“*  (3) 

2tt  ° 

where,  v » frequency 

Vq  = frequency  at  line  center 
Av  = full  width  of  the  line  at  half  maximum, 
and  the  relation 

/g(v)dv  * 1 (4) 

Is  satisfied.  Prom  equation  (3)»  It  Is  clear  that  the  value 

of  a quantity  as  a function  of  frequency,  P(v),  is  related  to 
Its  value  at  line  center,  , by 

ttAv 

P(v)  » — P^gCv)  (5) 

Prom  equations  (4)  and  (5),  one  obtains  the  expression 

itAv 

/P(v)dv  = P„  (6) 

2 ° 

which  relates  the  area  under  the  Lorentzlan  line  to  the  line- 
width  and  the  value  of  the  quantity  at  line  center. 

Lifetimes . An  atom  In  an  excited  state,  1,  can  sponta- 
neously decay  to  some  lower  state,  J,  by  emitting  a photon 
(radiative  decay)  or  a phonon  (nonradlatlve  decay).  If  the 
probability  of  a radiative  1-J  transition  Is  given  by  A^j , 
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then  the 


and  that  of  a nonradlatlve  1-J  transition  by  A"?^, 
probability  of  an  1-J  transition,  , is  given  by 


+A 


nr 

ij 


(7) 


where  all  the  transition  probabilities  are  expressed  in  sec 
The  lifetime  of  the  i-J  transition,  f , is  defined  as 


T 


ij 


(8) 


In  a similar  manner,  the  radiative  and  nonradlatlve  terms  in 
equation  (7)  can  be  expressed  in  terms  of  lifetimes  leading 
to  the  expression 


1 


1 1 
+ 

sp  ^ij 


(9) 


where,  » the  spontaneous  radiative  lifetime  of  the  i-J 
transition  = (A^j)~^, 

nr* 

^ij  " nonradlatlve  lifetime  of  the  1-J  transition 


^Aij;  . 


, th 


The  total  transition  probability  of  the  1 state  is  de- 
termined by  summing  equation  (7)  for  all  J.  The  lifetime  cor- 
responding to  this  total  transition  probability  will  be  refer- 
ed  to  here  as  the  fluorescent  lifetime,  Tf>  and  is  related  to 
the  transition  probability  by 


^f 


(10) 


Substituting  equation  (7)  into  (10)  yields 
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(11) 


— = Z CA^  +a"^) 

ZjLAij+AijJ 


which  can  be  written  as 


1 

T. 


1 1 


(12) 


'nr 


th 


state  and  t 


where  is  the  radiative  lifetime  of  the  1 . 

r nr 

Is  the  nonradiative  lifetime  of  the  1 state  corresponding 
to  (ZjA^j)“*  and  (ZjA^j)”'  respectively. 

Of  the  various  lifetimes  discussed  above,  only  x^  can 
be  measured  directly.  The  other  lifetimes  may  be  calculated 
from  the  relationships  presented  here  and  In  following  sec- 
tions . 

Cross  Sections . The  concept  of  cross  sections  is  often 
used  In  spectroscopy  as  an  alternate  means  of  expressing  tran- 
sition probabilities.  Some  of  the  more  Important  relation- 
ships involving  cross  sections  are  developed  below  using  the 
approach  found  In  DlBartolo  (Ref  24:406-411). 

The  absorption  cross  section  for  a transition  Is  defined 

by 

k(v) 


ct(v) 


(13) 


n, 


where  k(v)  » absorption  coefficient  at  frequency  v (cm”M, 
njj^  ■ population  density  of  the  lower  level  of  the 
transition  (cm~’). 

For  low  Intensity  excitations,  DlBartolo  (Ref  24:409)  has 
shown  that 
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I 


/k(v)dv  = (14) 

^sp 

where  X = transition  wavelength  In  the  material  (cm), 
g^  = degeneracy  of  the  upper  level, 

= degeneracy  of  the  lower  level, 
n^^  = population  density  of  the  lower  level  for  zero 
excitation  Intensity  (cm”’), 

Tgp  = spontaneous  radiative  lifetime  of  the  u-H  tran- 
sition (sec). 

Thus,  equation  (13)  can  be  rewritten  as 

X^  g 1 

/a(v)dv  — (15) 

^sp 

since  for  low  excitation  Intensities  n^^  = n®.  Using  the  re- 
sult of  equation  (6)  in  (15)  yields 

TT  X^  g 1 

— Ava  * (16) 

2 8’  'sp 

where  = a(v^).  Therefore,  the  peak  absorption  cross  sec- 


tion Is  given  by 


^ir’Av  g-  T 


I sp 


In  Nd:YAG,  the  degeneracies  of  the  various  Stark  levels  are 
the  same,  therefore,  the  Einstein  B coefficients  are  equal. 
Also,  since  Nd:YAG  Is  Isotropic,  T^p  can  be  replaced  by 
(Au^)”*  where  A^j^  Is  the  Einstein  A coefficient.  Substituting 


g„  - g-  and  T. 


(A^^^)”*  In  (17)  and  recalling  that 
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(18) 


yields 


B 


u£ 


ttXAv 


2h 


(19) 


Thus,  for  Nd:YAG,  the  absorption  cross  section  and  the  induced 
emission  cross  section  are  equal,  and  are  given  by  equation 
(17)  with  S^/g£  = 1. 

In  those  cases  where  the  absorption  coefficient  can  be 
measured  directly,  the  determination  of  a from  equation  (13) 
is  straight  forward.  However,  direct  absorption  measurements 
are  not  always  possible  and  an  alternate  method  of  determining 
cj  is  necessary.  One  such  method  depends  on  the  relationship 
between  the  induced  emission  cross  sections  and  fluorescent 
intensities  of  the  lines  in  the  fluorescence  spectrum  derived 
below. 

The  peak  emission  intensity  in  watts  of  the  i-j  transi- 
tion, > can  be  expressed  as 


(20) 


where,  * radiative  transition  probability  of  the  i-J 

transition, 

hVj^j  * energy  of  the  emitted  photon, 

= number  of  atoms  in  the  initial  state. 

Then  the  ratio  of  the  fluorescent  Intensities  of  i-J  and  1-k 
transitions  is  given  by 
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(21) 


Using  the  results  of  equation  (17)  with  the  ratio  of  the  de- 
generacies equal  to  one,  allows  equation  (21)  to  be  rewritten 
as 


fii. 

’■Ik  *"lk“lk''lk^ij  "l 


where  the  a's  are  the  peak  cross  sections.  If  the  Initial 
levels  are  in  thermal  equilibrium 


^1 

— = exp{-(E.-E.,  )/kT}  = a 

XT  ± ± 


(23) 


Substitution  of  equation  (23)  into  equation  (24)  and  some 
rearrangement  leads  to 


Ilk  ^ ^^ii^ii’^ii^ik^ii 
'ij  ^''ik'^lk'^ik^lj 


(24) 


where  n is  the  index  of  refraction.  Therefore,  if  the  emission 
cross  section  of  one  transition  is  known,  the  remaining  emis- 
sion cross  sections  can  be  determined  from  the  fluorescence 
data  and  equation  (24). 

Branching  Ratios . The  branching  ratio,  , is  a mea- 
sure of  the  fraction  of  the  total  fluorescence  emitted  into 
the  i-J  transition.  When  specified  in  terms  of  photon  rates, 
the  branching  ratio  is  given  by 


^i./^^l.1 


(25) 
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where  is  the  same  as  in  equation  (20).  If  one  takes  into 
account  the  finite  width  of  the  emission  line  and  substitutes 
for  the  frequency  in  terms  of  wavelength,  equation  (25)  becomes 

fXI. . (A)dA 

g = ^ (26) 

where  the  integral  in  the  numerator  is  over  the  i-J  emission 
line  and  the  term  in  the  denominator  is  the  sum  of  all  such 
Integrals.  Normally,  the  integrals  in  equation  (26)  are  de- 
termined from  the  fluorescence  spectrum  using  a planlmeter  or 
by  the  following  approximation.  Recall  that  for  Nd:YAG  the 
emission  lines  are  fairly  sharp  and  can  be  characterized  by 
the  Lorentzlan  llneshape  function.  For  sharp  lines  X is  es- 
sentially constant  over  the  line  and  equation  (26)  can  be  ap- 
proximated by 

X../I. .(X)dX 

g = y__il (27) 

where  X^j  is  the  wavelength  corresponding  to  the  peak  of  the 
emission  line.  For  a Lorentzlan  llneshape,  the  Integrals  in 
equation  (27)  can  be  evaluated  by  means  of  equation  (6)  in 
terms  of  peak  emission  intensities  and  linewidths.  Singh, 
et  al.,  have  used  both  of  the  above  methods  for  determining 
branching  ratios  in  Nd:YAG  and  obtained  excellent  agreement 
between  the  two  methods  (Ref  21). 

Another  important  relationship  involving  the  branching 
ratios  is  derived  below  after  the  method  of  Rushworth,  et  al. 
(Ref  25:3^-35).  From  equations  (20)  and  (25) 
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(28) 


For  a system  in  thermal  equilibrium,  N^,  the  population  of 
the  i^  state  is  related  to  N , the  population  of  the  ground 

O 

state  jy 


N 


j = Ng^xp(-(Ej-Eg)/kT)  = 


(29) 


Substituting  for  in  equation  (28)  yields 

^9 


ij 


, A j , a . 


(30) 


“i.j^lj^ig 

Since  the  '’P3/2  nianlfold  of  Nd:YAG  is  composed  of  two  Stark 
levels  that  are  in  thermal  equilibrium  (Ref  22:291))  the  sta- 
tistical average  of  the  total  radiative  transition  probability 
of  the  manifold  is  given  by 


A = 

^i^ig 


(31) 


Prom  the  relationship  between  lifetimes,  the  radiative  life- 
time of  the  ‘‘Fj/j  manifold  can  be  expressed  as 


“ r f 


a 

i,j"lj^ig 


(32) 


Dividing  equation  (30)  by  equation  (32)  yields 


or, 


ij 


■^r  ^i^lg 


(33) 


(34) 


T„a.  _ 

r Ig 


sp 
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Substitution  of  equation  (3^)  for  t in  equation  (17)  with 

sp 

gu/g£  = 1 gives 


^ ^ ^ii^  ^i^ig 

a ^TT^Ava. 

0 ig 


(35) 


where  is  the  peak  cross  section  of  the  i-J  transition  and 
Av  the  corresponding  linewidth.  Thus,  the  radiative  lifetime 
of  the  '*'^3^2  manifold  can  be  calculated  from  experimentally 
determined  values  of  the  branching  ratios,  linewldths,  and  cross 
sections . 

Quantum  Efficiency . The  term  quantum  efficiency  is  used 

differently  by  the  spectroscopist  and  the  laser  scientist.  In 

spectroscopy,  the  quantum  efficiency,  ri^,  is  usually  defined  by 

number  of  photons  emitted  from  i^  level 

rip  = (36) 

number  of  photons  absorbed  by  the  material 

and  is  a measure  of  the  overall  fluorescent  efficiency  of  a 

particular  level.  As  stated  previously,  in  Nd:YAG  each  pump 

photon  absorbed  relaxes  nonradlatively  to  the  ‘’F3/2  manifold, 

therefore,  can  be  expressed  as  the  ratio  of  the  radiative 

transition  probability  of  the  manifold  and  the  total 

transition  probability  of  the  manifold.  Therefore,  in  terras 

of  lifetimes 


n 


^f 


(37) 


When  the  laser  scientist  uses  the  term  quantum  efficiency 
he  is  most  often  referring  to  the  fluorescent  efficiency  of 
the  laser  transition.  That  is: 
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number  of  photons  emitted  in  laser  transition 


C38) 


number  of  photons  absorbed  by  the  material 
For  the  Nd:YAG  “Fj^j  manifold,  the  same  argument  used  In  ar- 
riving at  equation  (37)  yields 


Tf 


(39) 


"L  "R 

sp  sp 

The  values  of  and  ^.re  normally  calculated  from  ex- 
perimental data,  however,  direct  measurements  of  rij^  for  the 
manifold  of  Nd:YAG  have  been  made  (Refs  21,26). 


Spectroscopic  Parameters  of  Nd; YAG 

The  values  of  selected  spectroscopic  parameters  of  NdrYAG 

are  presented  in  Table  I.  The  values  of  and  are  for  the 

'*F,,,  manifold  and  those  of  and  t _ are  for  the  laser  tran- 

0 sp 

sltion.  The  multiple  values  of  a , t and  Po  of  Neeland  and 

0 sp  R 

Evtuhov  were  determined  by  different  methods  using  the  same 
sample.  Also,  Singh,  et  al.,  and  Kushida,  et  al.,  performed 
their  measurements  on  the  same  sample. 

It  is  clear  from  Table  I that,  with  the  exception  of  the 
fluorescent  lifetime,  there  is  little  agreement  among  the  var- 
ious investigators.  This  disagreement  has  created  a contro- 
versy over  the  role  of  nonradlatlve  processes  in  the  decay  of 
the  manifold  and  led  to  speculation  that  It  may  be  pos- 

sible to  affect  significant  improvements  In  the  overall  ef- 
ficiency of  Nd:YAG  as  a laser  material  (Ref  31).  This  contro- 
versy is  discussed  briefly  In  the  final  section  of  this  report 
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Table  1 

Spectroscopic  Parameters  of  Nd:YAG 


IV.  Power  Output  from  Nd : YAG  Lasers 


Tii^  laser  literature  is  full  of  expressions  for  predic- 
ting the  power  output  of  various  laser  systems  in  the  differ- 
ent operating  modes.  However,  In  the  majority  of  these  ex- 
pressions the  fundajnental  parameters  of  the  host  material  do 
not  appear  explicitly,  making  it  difficult  to  correctly  deduce 
the  functional  dependence  of  power  output  on  the  host  para- 
meters. In  this  chapter,  the  basic  expressions  for  a four 
level  laser  system  are  presented  in  terms  of  the  fundamental 
host  parameters  and  the  host  parameters  that  are  most  likely 
responsible  for  power  output  variations  in  Nd:YAG  are  identi- 
fied and  discussed. 

Nd : YAG  as  a Four  Level  Laser 

The  distinguishing  feature  of  four  level  lasers  is  the 
position  of  the  terminal  laser  level.  In  four  level  lasers, 
this  level  is  far  enough  above  the  ground  level  that  its  ther- 
mal population  with  respect  to  that  of  the  ground  level  is 
negligible.  Recall,  from  Fig.  3t) , that  in  Nd:YAG  the  termi- 
nal level  is  approximately  2000  cm~^  above  the  ground  state, 
therefore,  the  thermal  population  of  this  level  is  smaller 
than  that  of  the  ground  state  by  a factor  of  10“ ®.  Thus, 
Nd:YAG  can  be  modeled  using  a four  level  scheme  such  as  that 
Illustrated  in  Fig.  4.  In  this  somewhat  simplified  model,  R 
atoms  per  second  are  excited  from  the  ground  level  (1)  to  the 
pump  level  (4).  These  atoms  (in  the  case  of  Nd)  relax  non- 
radlatlvely  to  the  upper  laser  level  (3)  in  a time,  t^,  that 
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Pig.  4.  Schematic  of  a Four  Level 
Laser  Medium 


Is  short  compared  to  the  lifetime  of  the  upper  laser  level, 

Tf.  Relaxation  from  level  3 may  occur  by  a radiative  3-2 
transition  (laser  transition),  a nonradiatlve  3-2  transition, 
or  radiative  and  nonradiatlve  3-1  transitions.  Therefore, 
only  a fraction  of  the  atoms  in  level  3 contribute  to  the  3-2 
radiative  transitions.  This  fraction  corresponds  to  the  quan- 
tum efficiency,  In  the  case  of  Nd,  the  atoms  in  level  2 

undergo  nonradiatlve  transitions  to  the  ground  level  in  a 
time,  Tj,  that  is  short  compared  to  the  lifetime  of  level  3- 
Since  T2<<Tf«,  the  population  of  the  lower  level  is  essen- 
tially unchanged  from  its  thermal  equilibrium  value.  Thus, 
the  population  inversion,  AN,  varies  as  the  population  of 
the  upper  laser  level,  Nj . 

The  Gain  Coefficient 

Recall  that  induced  transitions  differ  from  spontaneous 
transitions  in  several  important  ways.  The  most  obvious  dif- 
ference is  that  induced  transitions  occur  in  both  directions 
while  spontaneous  transitions  are  always  emission  processes. 


Also,  spontaneous  emission  is  independent  of  the  field  inten- 
sity while  the  induced  processes  are  not.  This  dependence 
on  field  Intensity  is  reflected  in  the  induced  transition 
probability  rate,  W^(v),  which  from  Ref  32  and  Ref  33:82  is 
given  by 

W.  (v)  = ^(v)  (40) 

SirhvT^^ 

sp 

where,  = intensity  of  the  incident  field  at  frequency  v 

(watts/cm^ ) 

g(v)  = normalized  llneshape 

T = radiative  lifetime  of  the  transition, 
sp 

Plnlaily,  spontaneous  emission  produces  photons  that  have  no 
definite  phase  or  directional  relationship  with  the  incident 
field  while  induced  emission  produces  photons  that  are  coher- 
ent with  respect  to  the  incident  field  and  are  emitted  in  the 
same  direction  as  the  incident  field.  Thus,  each  photon  re- 
sulting from  Induced  emission  contributes  to  the  intensity  of 
the  incident  field  while  only  a small  number  of  the  spontaneous 
photons  do. 

Next,  in  the  manner  of  Ref  33:83,  consider  the  case  where 
a monochromatic  plane  wave  of  frequency  v and  intensity  is 
incident  on  a medium  having  n^  atoms  per  unit  volume  in  level 
3 and  n^  atoms  per  unit  volume  in  level  2.  Prom  the  preceding 
discussion  it  is  clear  that  the  only  significant  interactions 
between  the  medium  and  the  incident  field  are  the  result  of 
Induced  transitions.  Therefore,  the  power  per  unit  volume 
exchanged  by  the  field  and  the  medium  is  given  by 
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p 


(^1) 


= Cnj-Hj )W^Cv)hv 

Volume 

where  W^(v)  is  defined  in  equation  (40).  As  discussed  above, 
induced  transitions  are  coherent  with  respect  to  the  incident 
field,  so  equation  (4l)  may  be  interpreted  as  the  change  in 
intensity  per  unit  length  (dl^/dz).  Substituting  equation 


(40)  into  (4l)  yields 

dz 


(nj-nj ) 


X^g(v) 


Birr. 


sp 


(42) 


Equation  (42)  has  the  solution 


I^(z)  = I^(o )exp{Y (v)z} 


(43) 


where  I^(o)  is  the  initial  Intensity  and  y(v)  is  the  gain 
coefficient  defined  by 


y(v) 


(nj-n, ) 


X^g(v) 


Stti 


sp 


(44) 


Prom  equations  (43)  and  (44),  it  is  clear  that  when  nj  > 
the  incident  wave  is  amplified  and  laser  action  may  occur. 


Gain  Saturation  in  a Homogeneously  Broadened  Laser  Medium 

Prom  equation  (44)  it  appears  that  once  a given  population 
inversion  is  achieved  the  gain  coefficient  is  fixed.  However, 
that  is  not  the  case.  When  a field  is  present  the  inversion 
tends  to  decrease  from  its  zero  field  value  because  the  pop- 
ulation of  the  upper  level  is  greater  than  that  of  the  lower 
level  but  the  Induced  transition  probability  rate,  Wj^(v),  is 
the  same.  The  exact  dependence  of  the  populations  on  field 
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Intensity  can  only  be  determined  from  a rate  equation  analy- 
sis such  as  the  one  that  follows.  The  method  used  here  is 
similar  to  that  of  Yariv  (Ref  33:90-93). 

Consider  a four  level  system  such  as  that  depicted  in 
Pig.  4 except  that  now  a monochromatic  plane  wave  of  intensity 
and  frequency  v corresponding  to  the  3-2  transition  is  in- 
cident on  the  medium.  Therefore,  it  is  necessary  to  Introduce 
the  Induced  transition  probability,  W^(v),  to  account  for  the 
interaction  of  the  Incident  wave  and  the  medium.  The  result- 
ing rate  equations  for  levels  3 and  2 assuming  homogeneous 
broadening  are : 

dnj  nj 

= R (n3-n2)W.  (v)  (45) 

dt  P Tf  ^ 


and 


dn^ 

dt 


N nz 

+ (nj-na )W^(v)  


(46) 


Under  steady  state  conditions  these  equations  can  be  rewrit- 
ten as : 


n 


3 


Rp-^n^Wi(v) 

{(1/Tf)+W^(v)} 


(47) 


and 


n3{  (1/Tgp)^-W^(v) } 

{(1/t2)+W^(v)} 


(48) 


Solving  equations  (47)  and  (48)  for  (n^-n^)  yields: 
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nj-nj 


(49) 


1+{T^+(1-(T^/Tgp) )T2}W^(v) 


Now,  Tj,/Tgp  is  always  less  than  or  equal  to  unity  and  for  this 
model  Tj  is  short  compared  to  Therefore,  equation  (49) 

can  be  rewritten  as 


nj-nj 


p f _ 0 

l+TfWi(v)  l+T^W^(v) 


(50) 


where  An^  = ^p^f  value  of  the  Inversion  in  the  absence 
of  a field.  Substitution  of  equation  (40)  into  equation  (50) 
yields 


03-02 


An. 


l+I  /I 
V s 


(51) 


where,  I is  the  saturation  intensity  defined  as 
s 


STThV 


I = 
s 


_S£ 


X^g(v)T^ 


(52) 


and  corresponds  to  the  value  of  for  which  the  population  in- 
version has  decreased  to  half  its  zero  field  value. 

From  equations  (44)  and  (51)  the  gain  coefficient  can  be 
expressed  as 


y(v) 


An^X*g(v)  Yg(v) 


(53) 


The  decrease  in  the  gain  coefficient  with  predicted  by 
equation  (53)  is  termed  gain  saturation  and  (as  is  discussed 
in  the  following  section)  is  responsible  for  the  steady  state 
behavior  of  laser  oscillators. 
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Laser  Oscillation 


Consider  a laser  resonator  such  as  that  shown  In  Fig.  5 
where  the  laser  medium  has  a gain  coefficient,  y{v),  and  a 


Pig.  5.  A Simple  Laser  Resonator 
loss  coefficient,  a.  An  optical  field  with  an  Initial  Inten- 
sity Iq  will  have  an  Intensity  I,  after  one  round  trip  In  the 
resonator,  given  by: 

I =■  IoR^R2exp{2(Y(v)-a)Jl}  (54) 

where  R^  and  R^  are  the  power  reflectivities  of  the  mirrors. 

For  oscillation  to  occur,  the  minimum  Intensity  after  the  round 
trip  must  be  equal  to  1^.  This  occurs  when: 

RjR2exp{2(Y(v)-a)il}  (55) 

Substituting  Rj  = r^  and  R^  * r^  Into  equation  (55)  and 
solving  for  y{\))  yields: 

1 

Yt(v)  * a - _ln(rjrj)  (56) 

i 

where  the  subscript  't'  Indicates  that  this  Is  the  minimum  or 
threshold  value  of  the  gain  coefficient  for  oscillation  to 
occur.  Similarly,  from  equations  (44)  and  (56),  the  thresh- 
old value  for  the  population  Inversion  Is  defined  by: 
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(57) 


Sttt  I 

n = Cn3-n^l  = TTlV"  " 

A^g(v)  I 

For  the  case  where  the  frequency  of  the  Incident  field  corres- 
ponds to  the  center  frequency,  of  the  3-2  transition, 
substitution  for  from  equation  (3)  yields 

4ir^AvT  1 

n^  = 2.{a ^ln(r  r )}  (58) 

Recall  from  equation  (51)  that  the  population  inversion 
decreases  as  the  radiation  intensity  increases  due  to  the  in- 
crease in  induced  transitions.  However,  when  the  inversion 
falls  below  the  threshold  value  of  equation  (58),  laser  oscil- 
lation ceases,  the  radiation  Intensity  is  essentially  zero, 
and  the  population  inversion  begins  to  increase  again.  When 
the  Inversion  reaches  threshold,  oscillation  starts  again  and 
the  inversion  begins  to  decrease.  Due  to  this  saturation  be- 
havior, the  value  of  the  steady  state  population  inversion  is 
the  same  as  the  threshold  value.  Therefore,  at  steady  state 
each  atom  in  the  upper  laser  level  in  excess  of  the  threshold 
population  results  in  a laser  transition  in  order  to  restore 
the  population  inversion  to  the  threshold  value. 

Power  Output  for  CW  Operation 

Prom  the  preceding  discussion,  it  is  clear  that  for  CW 
operation  the  population  inversion  density  must  be  maintained 
at  the  threshold  value  given  by  equation  (58).  Thus,  the  pop- 
ulation density  of  the  upper  laser  level  must  be  maintained 
at 
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(59) 


^3  = 


rij  + 


n. 


For  an  ideal  four  level  system  the  thermal  population  of  the 
lower  laser  level,  n^,  is  essentially  zero  and  can  be  neglect- 
ed with  respect  to  n^ . This  assumption  Is  not  a good  one  for 
Nd:YAG,  since  for  one  atom  percent  doping  the  thermal  popu- 
lation of  the  lower  laser  level  Is  roughly  10^®  atoms/cm^  which 
Is  of  the  same  order  as  the  value  for  n^  given  by  Yariv  (Ref 
33:1^9).  However,  the  main  point  of  this  section  Is  to  deter- 
mine the  variation  of  power  output  as  a function  of  the  host 
parameters.  Neglecting  nj  with  respect  to  n^  does  not  change 
the  functional  dependence  of  the  power  output  on  the  host  para- 
meters and  does  result  In  a somewhat  simpler  expression  for 
power  output,  therefore,  the  assumption  that  n^  = 0 will  be 
made . 


Under  the  further  assumption  that  each  pump  photon  ab- 
sorbed excites  an  Ion  Into  the  upper  laser  level,  the  mini- 
mum photon  pumping  rate  is 


min 


^L^f 


(60) 


where  the  quantum  efficiency,  accounts  for  the  fact  that 
only  a fraction  of  the  exclt-d  Ions  produce  radiation  at  the 
laser  frequency  and  V Is  the  volume  of  the  material.  For  pump- 
ing rates  (R)  greater  than  laser  oscillation  occurs  and 

each  excited  ion  produces  a laser  photon  of  frequency  v.  Thus, 
the  power  output  for  CW  operation  is  given  by 


hv(R  - 


«mln^ 


(61) 
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This  equation  can  be  rewritten  In  terms  of  host  parameters  by 

means  of  equations  (60),  (58),  and  (39)  yielding 

4ir^Vv^A.vn^Tf_{a-(i/)l)ln(rir2  ) } 

P = hv(R ^ ) (62) 

where  n Is  the  Index  of  refraction  and  the  relation  X = c/nv 
has  been  used. 

Energy  Output  for  Q-Swltch  Operation 

In  contrast  to  the  relatively  simple  developement  that 
led  to  equation  (62),  the  developement  of  a similar  relation- 
ship for  Q-switch  operation  is  rather  complex.  The  approach 
used  here  Is  after  that  of  Wagner  and  Lengyel  (Ref  34:2040-42) 
and  assumes,  in  addition  to  the  assumptions  of  the  previous 
section,  that: 

1) .  Replenishment  of  the  inversion  by  the  pump  during 

the  life  of  the  pulse  is  negligible. 

2) .  No  change  in  the  inversion  occurs  during  the  switch- 

ing process  - l.e.  fast  switching. 

3) .  The  effects  of  continuous  pumping  and  spontaneous 

emission  on  the  inversion  are  neglected. 

Prom  Pig.  5,  the  optical  path  length  in  the  oscillator 
is  nt  where  n is  the  index  of  refraction  of  the  material  and 
the  time  of  transit  is 


t*  = nl/c  (63) 

Recall  that  the  cavity  lifetime  is  defined  as 

t^  = nJl/cL  (64) 
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where  L is  the  total  single  pass  loss  given  by 


L = ctA-lnCpjr^  ), 


(65) 


If  yl  is  the  gain  per  pass  and  (p  is  the  number  of  photons, 
then  the  rate  of  increase  in  photons  is  (pyi/t^  and  the  rate 


of  decrease  is  Thus 

d(ti 


— * {(yt-/t^)-(l/t  )}<t) 
dt  ° 


(66) 


For  each  photon  produced  the  total  inversion  population,  N, 
is  decreased  by  2,  therefore 


dN 

— = -2(rZ/t^)<t> 
dt 

Let  T = then  equations  (66)  and  (67)  become 

d(j) 

— = { (y/y^  )-1}4) 
dx 


(67) 


(68) 


and 


dN  Y 

— = -2—4. 
dx  Y,. 


(69) 


where  y^  = t^/Zt^  is  the  threshold  gain.  Since  y * N,  y/y^. 
can  be  replaced  by  N/N^  with  the  results 


and 


d4) 

— - {(N/N.  )-l}4. 
dx  ^ 


dN 

— = -2(N/N.  )(|) 
dx  ^ 


(70) 


Dividing  equation  (70)  by  (71)  and  integrating  yields 
<p  = ^{N^ln(N/N^)-(N-Nj^)} 
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(71) 


(72) 


*6 


since  the  initial  number  of  photons,  is  zero  if  sponta- 

neous emission  is  neglected.  Also,  substitution  of  equation 
(70)  into  (71)  and  some  rearranging  results  in 


d(N+2(|)) 


= -2^ 


The  number  of  photons  emitted  in  a single  pulse  is  obtained 
from  (73)  hy  integrating  over  the  duration  of  the  pulse  and 


is  given  by 


= /(f.dT  = ^{N^-N^+2((J)^-(|.^)} 


of  since  = 


N,-Nf 


An  expression  for  (N^-N^)  in  terms  of  the  threshold  inversion 
is  obtained  from  equation  (72)  by  letting  (p  = (p^  = 0 and 
N = with  the  result 

N^-N^  = N^ln(N^/N^)  (76) 

Now,  in  general  the  energy  stored  in  the  pulse  is 

N.  -Nf. 

E = — — (77) 

0 2 

where  v is  the  laser  transition  frequency.  However,  two  cases 
are  of  primary  Interest . 

Case  For  values  of  N^/N^  > 4,  Wagner  and  Lengyel 
(Ref  34:2042)  have  shown  that  >>  N^.  Therefore,  equation 
(77)  reduces  to 


= N.hv/2 
0 1 
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Case  II.  For  N /N  = 1,  N /N  = 1 (Ref  34;2042)  and 

XU  XX 


equation  (76)  can  be  written  as 


Nj,  ln(N^/N^)  ln(x) 
(N^/N^)-l  x-1 


(79) 


But 


ln(x)  = x-1  - — (x-l)^+...  for  2 > x > 0 

2 


Therefore 


or 


1 + 


N. 


2N. 


Ni  1 

— (N.-N.)  = -(N.-N.) 

N,  2 

u 

Substituting  equation  (82)  into  (77)  yields 


(80) 


(81) 


(82) 


Eq  = {(N|/N^)-N^}hv 


Since  = n^V/n^,  where  n^  is  given  by  equation  (58) 

N|X^Tf. 


= {- 


(a-(l/il)ln(rir2 ) ) 
sp 


- N. }hv 


or  in  terms  of  the  laser  transition  frequency 

- Nj^}hv 


Nfc^Tf. 

E = { i 

4Tr^v^n^AvTgp(a-(l/ll)ln(r  jr2 ) ) 


(83) 


(84) 


(85) 


At  first  glance,  it  might  seem  from  equation  (78)  that 
the  energy  in  the  pulse  is  Independent  of  the  host  parameters, 
however,  the  validity  of  the  assumption  Nj|_/N^  > 4 depends 
on  the  host  parameters  through  N^.  Thus,  the  energy  output 
is  indirectly  dependent  on  the  same  host  parameters  that  ap- 
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pear  explicitly  In  equation  (85)- 


Key  Host  Parameters 

The  host  parameters  that  may  affect  the  performance  of 
Nd:YAG  laser  systems  are,  from  Inspection  of  equations  (62) 
and  (85),  the  Index  of  refraction  (n) , the  loss  coefficient 
for  the  laser  transition  (ot),  the  llnewldth  (Av),  the  radia- 
tive lifetime  of  the  laser  transition  (t  ),  and  the  lifetime 

sp 

of  the  **P3/2  level  While  the  functional  dependence  of 

the  equations  on  these  parameters  Is  slightly  different.  It 
Is  clear  that  In  both  cases  Increases  In  n,  Av,  a,  and  t 
and  decreases  In  result  In  degradation  of  laser  performance. 

The  variations  of  most  of  the  host  parameters  with  tem- 
perature and  concentration  are  known.  However,  the  variations 
in  power  output,  of  Interest  here,  are  observed  between  rods 
with  nominally  the  same  concentration  and  operated  under  the 
same  thermal  conditions.  Therefore,  the  parameters  that  are 
mcie  sensitive  to  variations  in  concentration  and  temperature 
are  the  ones  of  primary  Interest. 

In  the  case  of  Nd:YAG,  the  Index  of  refraction  is  weakly 
dependent  on  concentration  for  normal  (=one  atom  percent) 
doping  levels  (Ref  12:91-92)  and  Is  essentially  constant  for 
small  changes  In  temperature  (Ref  8).  The  llnewldth  is  also 
a function  of  temperature  since  Nd:YAG  Is  thermally  broaden- 
ed. However,  for  small  changes  in  temperature  variations  In 
llnewldth  are  insignificant.  There  Is  no  evidence  to  Indi- 
cate that  the  llnewldth  Is  concentration  dependent  for  the 
normal  concentrations  used  in  Nd:YAG.  The  radiative  lifetime 
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I of  the  laser  transition  is  Independent  of  temperature  (Ref  35: 

1100)  and,  since  t is  determined  by  the  Interaction  of  the 

sp 

! fields  of  the  ion  and  the  crystal,  the  dependence  on  concen- 

I 

tratlon  should  be  negligible  except  at  extremely  high  doping 
levels.  The  loss  coefficient  varies  directly  with  the  absorp- 
tion coefficient,  k(v),  which  is  strongly  dependent  on  the 
population  of  the  terminal  laser  level  (recall  equation  (IM)). 
The  terminal  level  population  is  directly  proportional  to  the 
concentration  of  neodymium  ions  and  dependent  on  temperature 
through  the  Boltzmann  factor  (exp{-AE/kT} ) . Also,  the  loss 
coefficient  exhibits  strong  dependence  on  the  concentration 
of  impurity  ions  that  can  absorb  laser  radiation  and  the  amount 
of  scattering  in  the  crystal.  The  lifetime  of  the  upper  laser 
level  is  essentially  independent  of  temperature  (Ref  36:498), 
but  exhibits  a rapid  decrease  with  concentration  for  doping 
levels  above  one  atom  percent  (Ref  9 and  Ref  13:53). 

I Based  on  the  preceding  discussion,  the  key  host  parame- 

' ters  for  Nd:YAG  appear  to  be  a and  t^.  Indeed,  the  dlscus- 

» slon  of  energy  transfer  mechanisms  in  the  next  chapter  empha- 

i 

sizes  mechanisms  that  would  affect  these  parameters.  However, 

.1  experimental  correlations  between  the  variations  of  the  host 

parameters  and  laser  performance  have  not  been  made;  there- 
fore, the  possibility  that  another  parameter (s)  is  responsible 
for  the  variations  in  laser  performance  from  rod  to  rod  can- 
I not  be  neglected. 
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V.  Energy  Transfer  Mechanisms 


In  the  previous  chapter,  variations  in  laser  performance 
were  related  to  variations  in  the  parameters  of  the  host  mat- 
erial. Before  the  variations  in  these  parameters  can  be  con- 
trolled, the  mechanlsm(s)  that  causes  them  must  be  understood. 
The  discussion  that  follows  presents  qualitative  descriptions 
of  the  more  common  energy  transfer  mechanisms  and  their  ef- 
fects on  the  host  parameters,  as  well  as  comments  on  the  like- 
lihood that  a given  mechanism  is  responsible  for  variations 
in  these  parameters  from  one  rod  to  the  next. 

Resonant  Transfer 

Consider  the  case  where  the  metastable  “Fa/a  level  of  the 
neodymium  ion  coincides  with  an  energy  level  (above  the  ground 
state)  of  an  impurity.  These  levels  are  said  to  be  in  re- 
sonance and  energy  transfer  between  the  ions  may  occur.  Thus, 
a neodymium  ion  in  the  “Pa/a  level  may  transfer  its  energy  to 
the  impurity  ion  producing  a neodymium  ion  in  the  ground  state 
and  an  excited  impurity  ion.  As  illustrated  in  Pig.  6a,  the 
excited  impurity  ion  may  then  decay  either  radlatively  or  non- 
radlatlvely  to  its  ground  state. 

The  net  effect  of  resonant  transfer  is  to  deplete  the 
population  of  the  metastable  level  with  no  contribution  to  the 
fluorescent  yield  of  the  laser  transition.  Since  this  mech- 
anism represents  an  additional  decay  process,  it  should  re- 
sult in  a decrease  in  the  fluorescent  lifetime  of  the  ‘‘Pa /2 
level.  Furthermore,  if  the  impurity  ion  decays  by  radiative 
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Pig.  6.  Energy  Transfer  Mechanisms:  (a)  Resonant 
Transfer,  (b)  Radiation  Trapping,  (c)  Cross  Relaxa- 
tion, (d)  Auger  Recombination. 


means,  it  may  be  possible  to  identify  it  from  the  fluorescence 
spectrum  of  the  material.  The  strength  of  the  resonant  trans- 
fer process  depends  on  the  impurity  concentration.  It  is 
highly  unlikely  that  the  Impurity  concentration  is  the  same 
for  all  rods,  which  could  account  for  rod  to  rod  variations. 

Radiation  Trapping 

The  term  radiation  trapping,  as  used  here,  refers  to  the 
process  (see  Pig.  6b)  whereby  an  emitted  photon  encounters 
an  impurity  ion  with  a matching  transition  or  a lattice  de- 
fect and  is  reabsorbed. 

Since  radiation  trapping  does  not  involve  an  interaction 
with  the  metastable  level  of  the  neodymium  ion,  it  does  not 
affect  the  fluorescent  lifetime.  However,  this  mechanism  will 
Increase  the  loss  coefficient  of  the  material.  This  mechanism, 
like  resonant  transfer,  is  a likely  one  for  explaining  the 
performance  variations  among  rods  since  it  too  depends  on  the 
concentration  of  impurity  ions  or  defects.  Belt,  et  al. , 

(Ref  37)  have  found  that  substantial  amounts  of  iron  are  pre- 
sent in  the  insulating  materials  used  in  the  growth  process. 
Since  iron  exhibits  significant  absorption  at  l.Oby,  it  should 
be  a prime  suspect  as  the  cause  of  radiation  trapping  if  ab- 
sorption data  indicate  the  possibility  that  this  mechanism  is 
active  in  NdrYAG  (Ref  37). 

Cross  Relaxation 

Cross  relaxation  is  a process  somewhat  similar  to  rad- 
iation trapping  by  an  impurity  in  that  it  requires  matching 
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transitions.  However,  in  the  case  of  radiation  trapping,  the 
emitted  photon  is  rebsorbed,  while  for  cross  relaxation,  the 
Interaction  of  ions  with  such  matching  transitions  results  in 
an  entirely  different  relaxation  process.  As  an  example  of 
cross  relaxation,  consider  the  case  (see  Pig.  6c)  where  an 
excited  neodymium  ion  interacts  with  another  neodymium  ion  in 
the  ground  state.  The  excited  ion  may  undergo  a radiative 
transition  to  one  of  the  ‘*1  manifolds  above  the  state 

followed  by  a nonradiative  transition  to  the  ground  state, 
while  the  unexcited  ion  undergoes  a corresponding  absorption 
from  the  ground  state  and  subsequent  nonradiative  relaxation 
back  to  the  ground  state.  The  net  result,  is  depletion  of 
the  excited  state  population  with  no  contribution  to  the  fluor- 
escence at  the  laser  frequency.  The  probability  that  cross 
relaxation  will  occur  depends  on  the  proximity  of  the  inter- 
acting ions  (concentration)  as  well  as  the  energy  match  be- 
tween the  transitions.  In  those  cases  where  the  transitions 
do  not  exactly  match,  cross  relaxation,  with  the  absorption 
or  emission  of  a phonon  for  energy  conservation,  may  still 
occur  (Ref  24:497). 

Since  cross  relaxation  is  an  alternate  decay  mechanism 
for  the  metastable  state,  it  decreases  the  fluorescent  life- 
time. In  fact,  Danlelmeyer  has  proposed  that  phonon-assisted 
cross  relaxation  Involving  the  s/2  5/2 

transitions  is  responsible  for  the  decrease  in  the  fluorescent 
lifetime  of  the  '*P3/2  manifold  with  neodymium  concentration 
that  is  observed  in  Nd:YAG  (Ref  38).  Singh  (Ref  39)  has  shown 
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that  the  mechanism  for  cross  relaxation  proposed  by  Danielmeyer 
involves  the  absorption  of  a phonon  from  the  lattice  and  v;ould 
result  in  the  fluorescent  lifetime  exhibiting  strong  temper- 
ature dependence.  Since  the  fluorescent  lifetime  is  observed 
to  be  essentially  Independent  of  temperature  variations  (Ref 
36:498),  Singh  has  proposed  an  alternate  mechanism  for  cross 
relaxation  involving  the  3/2  5/2 

sltlons  and  the  emission  of  a phonon.  This  mechanism  predicts 
a decrease  in  fluorescent  lifetime  with  neodymium  concentra- 
tion, but  essentially  no  change  in  lifetime  with  temperature. 

Recall  that  concentration  gradients  of  20-30  percent  for 
20  cm  long  boules  are  known  to  exist  in  quality  boules  of 
Nd:YAG.  In  the  growth  of  Nd:YAG  at  concentrations  near  one 
atom  percent,  these  gradients  may  produce  Nd  concentrations 
high  enough  to  decrease  the  fluoresent  lifetime  of  the  meta- 
stable level  thru  cross  relaxation.  Therefore,  cross  relax- 
ation is  another  likely  cause  of  performance  variations  in 
Nd:YAG  rods. 

.Multiphonon  Nonradiative  Decay 

In  addition  to  radiative  decay  involving  photons,  excit- 
ed ions  may  relax  to  the  ground  state  by  the  nonradiative  em- 
ission of  phonons.  When  the  energy  gap  for  a transition  ex- 
ceeds that  of  the  highest  energy  phonons  in  the  material,  non- 
radiative decay  may  still  occur  by  the  emission  of  more  than 
one  phonon.  The  transition  probability  for  these  multiphonon 
emission  processes  is  given  by 
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Pi 


(86) 


w 

0 


exp (ftw^/kT ) 
.exp(ftaj^/kT)-lJ 


where  Is  the  transition  probability  for  T*0,  is  the 
phonon  energy,  and  p^  is  the  order  of  the  process  determined 

by 


P 


i 


AE  of  the  transition 


(87) 


It  has  been  shown  that  raultlphonon  decay  generally  involves 
phonons  that  result  in  the  lowest  order  process  consistent 
with  the  cut-off  frequency  of  the  phonon  spectrum  (Ref  40:35- 
55). 

If  multiphonon  decay  is  significant  in  a material,  the 
fluorescent  lifetime  should  decrease  with  increasing  tempera- 
ture. Recall  that  the  fluorescent  lifetime  of  the  metastable 
level  in  Nd:YAG  is  essentially  constant  with  temperature  (Ref 
36:498).  Therefore,  multiphonon  relaxation  is  probably  not  an 
important  decay  mechanism  in  Nd:YAG. 


Excited  State  Pumping 

If  a material  is  pumped  hard  enough,  the  population  of 
the  excited  state  will  eventually  reach  the  point  where,  if  a 
suitable  transition  exists,  the  pumping  radiation  is  nearly  as 
likely  to  be  absorbed  by  an  ion  in  the  excited  state  as  by  an 
ion  in  the  ground  state.  This  excited  state  absorption  dep- 
letes the  population  of  the  metastable  level  and  makes  it  more 
difficult  to  achieve  or  maintain  a population  Inversion.  If 
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I. 


(as  in  Nd:YAG)  the  levels  above  the  inetastable  level  relax 
nonradlatively  back  to  the  metastable  level,  excited  state 
pumping  does  not  deplete  the  Inversion.  However,  the  effici- 
ency of  the  pump  In  populating  the  metastable  level  Is  decreas- 
ed. Continued  Increases  in  pump  power  result  in  further  de- 
creases in  pump  efficiency  and  saturation  of  the  fluorescence 
from  the  metastable  level. 

Since  the  energy  of  the  metastable  level  of  Nd  is  much 
greater  than  that  of  the  ground  level,  the  population  of  the 
metastable  level  is  much  less  than  that  of  the  ground  level 
even  during  laser  operation.  Therefore,  excited  state  pumping 
is  probably  not  significant  in  Nd:YAG. 

Auger  Recombination 

If  two  ions  in  the  metastable  level  interact,  it  is  pos- 
sible that  one  ion  will  lose  its  excitation  energy  to  the  oth- 
er producing  an  ion  in  a still  higher  state  and  an  ion  in  the 
ground  state.  If  the  doubly  excited  ion  relaxes  nonradlative- 
ly to  the  metastable  level,  the  net  effect  of  this  process  is 
to  remove  one  ion  from  the  excited  state  population  with  no 
contribution  to  the  fluorescent  output.  This  process,  illus- 
trated in  Pig.  6d,  is  called  Auger  recombination  (Ref  38). 

Auger  recombination,  is  an  effective  fluorescence  quench- 
ing process  only  for  pump  powers  beyond  those  normally  used 
in  laser  operation  (Ref  38).  Therefore,  Auger  recombination 
probably  does  not  play  a significant  role  in  laser  performance 
under  normal  operating  conditions. 


44 


VI.  Reconimendatlons  for 


Future  Investigations 

Prom  the  preceding  discussion.  It  should  be  clear  that, 
although  quite  a bit  Is  knovm  about  the  chemistry,  growth, 
and  spectroscopy  of  NdcYAG,  there  are  still  many  questions 
that  need  to  be  answered.  Some  future  investigations  that 
may  shed  light  on  the  mechanism  or  mechanisms  responsible  for 
the  output  power  variations  in  NdcYAG  are  discussed  in  the 
following  pages . 

At  the  conclusion  of  Chapter  III,  it  was  mentioned  that 
the  discrepancies  in  the  values  of  the  spectroscopic  parame- 
ters of  NdcYAG  had  led  to  speculation  that  significant  im- 
provements in  the  laser  performance  of  NdcYAG  might  be  pos- 
sible. Recent  experiments  by  Singh,  et  al.  (Ref  39),  in  which 
single  phase  powders  of  NdcYAG  were  produced  which  had  signi- 
ficantly longer  fluorescent  lifetimes  than  commercial  NdcYAG 
of  the  same  concentration,  seem  to  support  Singh's  original 
assertion  that  the  correct  value  for  the  radiative  quantum 
efficiency,  in  commercial  NdcYAG  is  0.56  (Ref  21)  rather 
than  1.0  (Refs  22,26).  If  Singh's  value  for  in  commercial 
NdcYAG  is  correct,  then  it  may  be  possible  to  improve  the  ef- 
ficiency of  NdcYAG  by  a factor  of  two.  Now,  the  correct  val- 
ue of  whether  it  is  0.56  or  1.0,  has  no  bearing  on  rod  to 
rod  performance  variations  if  it  is  constant.  However,  the 
value  of  does  provide  important  information  on  the  role  of 
nonradlatlve  processes  in  the  decay  of  the  level  which 
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may  help  In  determining  the  mechanism  or  mechanisms  responsible 
for  performance  variations  in  Nd:YAG.  Therefore,  one  goal  of 
future  investigations  should  be  to  resolve  the  controversy 
over  the  value  of  in  commercial  Nd:YAG. 

If  the  mechanlsmCs)  responsible  for  output  variations 
in  Nd:YAG  is  ever  to  be  determined,  then  experimental  data 
relating  laser  performance  and  variations  in  the  host  para- 
meters are  essential.  Therefore,  one  of  the  immediate  goals 
of  any  further  investigations  in  this  area  should  be  to  pro- 
vide that  data.  One  possible  approach  for  doing  this  is  out- 
lined below. 

Any  investigation  into  the  relationship  between  host  pa- 
rameters and  laser  performance  must  use  rods  for  which  the 
complete  growth  history,  including  the  position  of  the  rod  in 
the  boule,  is  known.  The  rods  employed  should  have  concentra- 
tions of  approximately  0. 8-1.0  atom  percent  (the  normal  range 
for  Nd’.YAG).  All  rods  should  have  the  same  dimensions  and 
external  finish.  The  CW  laser  output  of  each  rod  should  be 
determined  using  a standard  laser  cavity  and  Identical  pump- 
ing conditions.  Once  the  laser  performance  of  the  rods  has 
been  determined,  the  following  passive  measurements  should  be 
made : 

1) .  Determine  the  relative  optical  quality  of  the  rods 

using  Twyman-Green  interferometry. 

2) .  Measure  the  passive  absorption  from  0. 3-2.0  microns 

at  room  temperature  and  at  approximately  the  same 
temperature  as  that  experienced  by  the  rod  during 
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laser  operation. 

3) .  Measure  as  a function  of  distance  along  the  rod. 

I 

4) .  Obtain  the  fluorescence  spectrum  of  the  '*^'’3/2 

fold  and  determine  the  branching  ratios  and  laser 
transition  cross  section. 

5) .  Measure  directly  (Refs  21,26). 

6) .  Determine  crystal  defect  concentration. 

7) .  Determine  the  actual  neodymium  concentration  as  a 

function  of  distance  along  the  rod. 

8) .  Analyze  the  rods  for  impurities  such  as  iron  that 

absorb  the  laser  or  pump  wavelengths. 

Clearly,  an  Investigation  of  the  scope  of  the  one  Just  out- 
lined is  expensive  and  time  consuming  and  may  never  be  prac- 
tical from  an  economic  standpoint.  However,  unless  the  num- 
ber of  possible  mechanisms  can  be  narrowed  down  substantially, 
the  above  approach  may  be  necessary. 

Based  on  the  available  data,  one  possible  explanation 
for  the  rod  to  rod  variations  seems  to  deserve  special  con- 
sideration. Recall  that  the  low  distribution  coefficient  of 
neodymium  leads  to  melt  enrichment  and  produces  concentration 
gradients  in  the  boule.  Since  the  fluorescent  lifetime  is 
known  to  decrease  rapidly  for  concentrations  above  1.2  atom 
percent,  concentration  gradients  can  result  in  decreased  laser 
performance.  Furthermore,  the  reported  values  of  the  distri- 
bution coefficient  range  from  0.12-0.21,  thus,  the  same  start- 
ing melt  composition  and  growth  conditions  may  not  always  pro- 
duce the  same  concentration  profile  in  the  boule.  The  net 
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result  is  that  it  is  possible  that  a boule  grown  to  contain 
a specific  neodymium  concentration  may,  in  fact,  contain  a 
significantly  higher  concentration  than  desired  and  have  a cor- 
respondingly lower  fluorescent  lifetime  with  decreased  laser 
performance . 

Experimental  confirmation  of  the  preceding  hypothesis 
would  involve  the  following: 

1) .  Boules  should  be  grown  from  the  same  melt  composi- 

tion and  under  the  same  growth  conditions. 

2) .  Rods  should  come  from  the  top,  middle,  and  bottom 

of  the  boules. 

3) .  All  rods  should  have  the  same  dimensions  and  exter- 

nal finish. 

4) .  Evaluate  laser  performance  based  on  CW  power  out- 

put in  a standard  cavity. 

5) .  Compare  optical  quality  of  rods  using  Twyman-Green 

Interferometry . 

6) .  Measure  as  a function  of  distance  along  the  rod. 

7) .  Determine  neodymium  concentration  as  a function  of 

distance  along  the  rod. 

If  the  hypothesis  is  correct,  the  rods  which  exhibit  the  poorer 
performance  will  have  shorter  fluorescent  lifetimes  and  higher 
concentrations . 
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